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Abstract: A campaign of 35 epochs of milli-arcsecond resolution VLBA observations of the archetype colliding - 
^ . wind WR+O star binary system WR 140 show the wind-collision region (WCR) as a bow-shaped arc of emis- 
^| sion that rotates as the highly eccentric orbit progresses. The observations comprise 21 epochs from the 1993- 
2001 orbit, discussed by Dougherty et al. (2005), and 14 epochs from the 2001-2009 orbit, and span orbital 
phase 0.43 to 0.95. Assuming the WCR is symmetric about the line-of-centres of the two stars and "points" at 
the WR star, this rotation shows the O star moving from SE to E of the WR star between these orbital phases. 
Using IR interferometry observations from IOTA that resolve both stellar components at phase 0.297 in con- 
junction with orbital parameters derived from radial velocity variations, the VLBA observations constrain the 
inclination of the orbit plane as 120° ± 4°, the longitude of the ascending node as 352° ± 2°, and the orbit semi- 
major axis as 9.0±0.1 mas. This leads to a distance estimate to WR 140 of 1.81 ±0.08 kpc. Further refinements 
O ■ of the orbit and distance await more IR interferometric observations of the stellar components directly. 



1 Introduction 

The 7.9-year period WR+O system WR 140 (HD 193793) is the archetype of CWB systems. It is 
comprised of a WC7 star and an 04-5 star in a highly elliptical orbit (e ~ 0.88), where the stellar 
separation varies between ~ 2 AU at periastron to ~ 30 AU at apastron. This highly eccentric 
orbit clearly modulates the dramatic variations observed in the emission from the system, from X- 
ray energies to radio wavelengths (Williams et al. 1990). Perhaps the most dramatic variations are 
observed at radio wavelengths, where there is a slow rise from a low state close to periastron of a 
few mly, to a frequency-dependent peak in emission of 10's of mly between orbital phase 0.65 to 
0.85, before a precipitous decline immediately prior to periastron (see Fig 1). A number of attempts 
to model these variations have been made (e.g. Williams et al. 1990, White & Becker, 1995) with 
limited success, though advances in our understanding of WCRs are being made (e.g. Pittard & 
Dougherty (2006); Pittard, these proceedings). Accurate orbital parameters are critical inputs to these 
models. 

Many of the orbital parameters in WR 140, in particular the orbital period (F), epoch of periastron 
passage (T D ), eccentricity (e), and the argument of periastron (u) have been established by others (see 
Marchenko et al. 2003 and references therein), and refined most recently in an extensive observing 
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Figure 1: Radio emission from WR 140 at 15 GHz (blue circles), 5 GHz (green triangles), and 
1.6 GHz (red squares) as measured with the VLA (solid) and WSRT (open). Data are from Williams 
et al. (1990) and White & Becker (1995). 

campaign during the 2009 periastron passage. However, the orbital inclination (i), semi-major axis 
(a), and the longitude of the ascending node (ft) require the system to be resolved into a "visual" 
binary. The two stellar components in WR 140 have been resolved using the Infra-red Optical Tele- 
scope Array (IOTA) interferometer at a single epoch (Monnier et al. 2004). This single observation 
sets the scale and orientation of the orbit since it constrains potential families of solutions for (i, a, ft). 
Further epochs of IOTA observations have been completed, but until analysis is complete, the VLBA 
observations of the WCR offer the only means to determine the orbit direction and constrain i, and 
hence ft and a. 

An initial analysis of 21 epochs of VLBA observations taken between 1999 and 2000 (orbital 
phase 0.74 to 0.95) was described in Dougherty et al. (2005). The work presented in this paper is 
an amalgamation of those observations with an additional 14 observations obtained between 2004 
to 2008, that extended the orbital phase coverage from 0.43 to 0.96. A re-analysis of the earlier 
observations with the new data leads to tighter constraints on the derived orbit parameters. 

2 VLBA observations of WR 140 

WR 140 was observed with the VLBA at 8.4 GHz at 35 epochs between orbital phase 0.43 and 0.97. 
A selection of images are shown in Fig. 2. The 8.4-GHz emission is clearly resolved, with a bow- 
shaped emission region observed at many epochs. This shape is anticipated for a WCR from model 
calculations (e.g. see Eichler & Usov 1993, Dougherty et al. 2003, Pittard et al. 2006). The WCR 
rotates from "pointing" NNW to W over the observed orbital phases. This rotation is key to deter- 
mining the orbit of WR 140. Assuming the arc of emission is symmetric about the line-of-centres and 
points towards the WR star, the O star is to the SSE of the WCR at epoch 0.43, and approximately to 
the E at phase 0.96. 
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Figure 2: VLBA 8.6-GHz images of WR140 at phases 0.62, 0.74, 0.78, 0.87, 0.92, and 0.93 from the 
1993-2001 orbit (phases 0.74, 0.87 and 0.93 - taken from Dougherty et al. 2005) and the 2001-2009 
orbit. The synthesized beam is 2.0 x 1.5 mas 2 in the 1993-2001 orbit, and approximately 1.3x that 
for the latest observations. Note change in RA and Dec between the images taken from the 1993-2001 
orbit and the recent orbit. Rotation and proper motion of the WCR are evident during both orbits. 



Deriving the Orbit 

On June 17, 2003 Monnier et al. (2004) observed WR 140 to have a separation of 12. 91^4 mas at 
a position angle of 151. 7^3 degrees east of north. Using P = 2896.6 days, T Q = 2446156.3, 
e = 0.897 and u = 46.8° determined from analysis of observations during the 2009 periastron (Fahed 
et al. 2010), the observation at orbital phase 0.296 restricts potential sets of solutions for i, a and Q to 
those shown in Fig. 3 for inclinations in the range 0° < % < 180°. 

The change in the orientation of the WCR with orbital phase gives the inclination since each (i, ft) 
solution family provides a unique set of position angles as a function of orbital phase. A weighted 
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Figure 3: Solutions for the longitude of the ascending node (red line) and orbit semi-major axis (green 
line) as a function of orbit inclination, derived from an IOTA separation and a position angle of the 
stellar components at orbit phase 0.296 (Monnier et al. 2004). The uncertainty in the IOTA observation 
gives an error in ft of closely ±1°, and in the semi-major axis of ±0.3 mas. 

minimum x 2 measure between the observed position angle of the line of symmetry of the WCR, and 
by proxy the line-of-centres of the stars, as a function of orbit phase determined for different sets of 
(i, ft) leads to a best-fit solution of % = 120° ± 4° and ft = 352° ± 2° (Fig. 4). These values lead to a 
semi-major axis of a = 8.97 ± 0.13 mas. 
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Figure 4: The change in the position angle of the axis of symmetry of the WCR as a function of orbital 
phase. The green line is the weighted best-fit curve of the position angle of the line-of-centres of the 
two stars as a function of phase, corresponding to i = 120° and ft = 352°. The other two lines show 
the fits for Ax 2 = ±1, showing the range of potential fits for the quoted uncertainties. 




Distance of WR 140 



Stellar distance is one of the more difficult parameters to determine for stars. Marchenko et al. (2003) 
determined a sin i = 14.10 ± 0.54 AU from radial velocity observations. From our estimate of i, this 
gives a = 16.28 ± 0.81 AU. A semi-major of a = 8.97 ± 0.13 mas then gives a distance estimate of 
1.81 ± 0.09 kpc, consistent with the previous estimate of 1.85 ± 0.16 kpc by Dougherty et al. (2005). 

Summary 

High-resolution radio observations of the WCR in WR140 have provided, to date, the only way to 
determine the direction and orientation of the orbit, starting from the scale of the orbit as deduced from 
IR interferometry and orbital parameters from optical spectroscopy. From this work it is possible to 
determine a precise, and hopefully accurate, distance to WR140. 

WR140 is the best colliding-wind binary system for attempts to understand the underlying par- 
ticle acceleration processes and physics in wind-collision regions, in large part due to the wealth of 
observational constraints. However, the orbit and distance are critical to modelling the WCR, and 
further refinements of these parameters await further IR interferometric observations that will resolve 
the stellar components directly. 
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